Level-spacing distribution of a fractal matrix 
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We diagonalize numerically a Fibonacci matrix with fractal Hilbert space structure of dimen- 
sion d/ = 1.8316... We show that the density of states is logarithmically normal while the corre- 
sponding level-statistics can be described as critical since the nearest-neighbor distribution function 
approaches the intermediate semi-Poisson curve. We find that the eigenvector amplitudes of this 
matrix are also critical lying between extended and localized. 

Pacs numbers: 71.30.+h, 61.44.-n, 71.23.Fb 

The numerical diagonalization of one-electron Hamiltonians in discrete space tight-binding lattices has been proved 
a very useful tool to treat the effect of disorder on the motion of quantum particles, alternatively to field theoretic 
methods. In particular, the direct computation of eigenvalues and eigenvectors combined with finite-size scaling 
techniques can answer questions about Anderson localization due to disorder [Q, quantum chaos in the energy levels 
or wave functions [^,^, etc. The difficulties become immense when one wants to treat many-body problems in the 
same framework since the dimension of the Hilbert space increases dramatically (exponentially) with the number of 
electrons. However, also in this case the nearest-neighbor hopping bounded kinetic energy and the two-body character 
of the interaction guarantee that the Hamiltonian matrix structure is very sparse, often described as multifractal in 
the adopted Hilbert space. For example, in ref. Q multifractal exponents Dq were computed to characterize the 
Hilbert space structure of interacting many-body Hamiltonians. 

We shall treat a much simpler, albeit interesting problem, with a given Hamiltonian which consists of zeroes and 
ones, neither periodic nor random but with a simple fractal structure. We construct and diagonalize a Fibonacci 
matrix of order n of size Nn x 7V„, with Nn — Nn-i + Nn-2, Nq = 7V_i = 1. This matrix represents a self-similar 
fractal object itself and was introduced in ref. jsj in order to obtain the ground state of the square Ising antiferromagnet 
in a maximum critical field. The Fibonacci matrix is a transfer matrix connecting the possible ground states of the 
n X m and n x (m + 1) square lattices. For size Nn it has the block form 



Fn-1 Gn-1 

G„ _ 1 



(1) 
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where Gn-i represents the iV„_i x iV„_2 submatrix of F„_i (its upper left corner), G'^_i its transpose and denotes 
the zero Nn-2 x Nn-2 matrix. The sequence of fractal matrices for n ~ 1, 2, ... begins with the matrix of size Ni — 2 
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etc and the corresponding matrices for higher n's are easily obtained from eq. (|i|) 
matrix Fiq of size iVio — 144 is shown in fig. 1. 



successively. For example, the 
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FIG. 1. The fractal Fibonacci matrix _Fio of size A^io = 144. The ones are denoted with black and the zeroes with white. 



In the Fibonacci matrix structure of fig. 1 we observe no coverage of the full two-dimensional Hilbert space, with 
blocks within blocks, so that the matrix is a simple fractal object. The fractal dimension of the Fibonacci matrix can 
be obtained form eq. (1) if we associate the mass M„ measured by the number of I's to the Fn matrix of size iV„, via 
the scaling relation 

Mn=Nn'. (5) 

The corresponding asymptotic (n oo) size ratio Nn+i/Nn = (1 + \/5))/2 and mass ratio M„+i/M„ = 1 + \/2 for 
two successive generations n and n + 1 can enter eq. (5). The size ratio is obtained from two successive Fibonacci 
approximants (equals the ratio of the sides of two successive rectangles) and by definition corresponds to the golden 
mean (1 + \/5)/2. The mass ratio M„+i/M„ = 1 + \/2 is obtained from the of the series M„ = Af„_i + 2(M„_2 + 
M„_3 + Mn-i + ... + Afi + 2), n = 2, 3, .., Mi = 3 and the formation of Af„+i/M„ = 1 + 2( jg-yij— ^ + a/„/m„_, + 

•■■ + M /Afi ~^ aTT^) ^hich in the limit oi n ^ oo where = — > x, ^ = ^f^"^ m"~2 ~^ becomes 

X — \ + 2(i + ^ + ■■■) and by summing the geometric series we obtain x = 1 + \/2. Therefore, the matrix of fig. 1 is 
a simple fractal (monofractal) object with dimension between 1 and 2 given by 

log{l + V2) 



log{{l + y5)/2) 



1.8316... (6) 



In order to obtain the residual entropy of the antiferromagnetic Ising model in maximum critical field the leading 
eigenvalue of such Fibonacci matrices up to -Fg was estimated in ref. 

Our purpose is entirely different. We shall consider the full spectrum of the Fibonacci matrix (eigenvalues and 
eigenvectors) and address questions related to Anderson localization and quantum chaos. We emphasise that we do 
not have in mind any particular physical system described by this matrix although such fractal matrices might be 
in general relevant for critical (between extended and localized) one-electron states in disordered or quasiperiodic 
potentials or possibly many-electrons states for interacting electrons in the presence of disorder . We shall address 
the following questions concerning the Fibonacci matrix: 

(1) What is the density of states? 

(2) What is the required 'unfolding' procedure and the nearest-level spacing distribution function? 
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(3) What is the behavior of the corresponding eigenvectors? 
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FIG. 2. (a) The density of states (full circles joined by a solid line) for the Fibonacci matrix Fiq of size A'^ = A^ig = 10946. (b). The 
corresponding density of states (full circles joined by a solid line) as a function of log{E), where E is the absolute value of the energy. The 
simple Gaussian fit with mean —1.79 and standard deviation 0.91 is also shown by the continuous line. 



Firstly, we compute the density of states p{E) as a function of the energy E which is shown in fig. 2(a). It is seen to 
be dominated by a strong peak at the origin E = and also has very long tails which increase with the size. Instead, 
the distribution of p{logE) shown in fig. 2(b), where E is the absolute value of the energy, leads to a Gaussian-like 
function. From the Gaussian fit of fig. 2(b) we can conclude that the density of states p{E) of the Fibonacci matrix 
follows a log-normal distribution. We have obtained similar densities of states for all other sizes examined. Having 
answered the first question, we proceed to compute the energy level-statistics in this system. In order to do so it is 
required to solve the problem of 'unfolding' for the unusual density of states of fig. 2 which demonstrates a strong peak 
at the origin. It is impossible to obtain a meaningful level-statistics without solving the 'unfolding' problem which is 
required to make the density of states (particularly this unusually singular density) constant, with mean- level spacing 
A = 1. Our approach consists in 'unfolding' the energy levels by using the density of fig. 2(b) instead of using that 
of fig. 2(a). In fact, this is essentially a double unfolding procedure of the original spectrum, one going from fig. 2(a) 
to 2(b) and the second done for the data of fig. 2(b), so finally the mean level-spacing becomes A = 1. 
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FIG. 3. (a), (b) The calculated nearest level-spacing distribution P{S) by unfolding the data for p(log{E)) for two Fibonacci matrices 
of sizes N = Ni7 = 4181, A^ig = 10946. The dashed Une is the Wigner surmise P{S) = (tt / 2) Sexp{- (17/4)3^) (extended states), the 
dotted line is the Poisson law P{S) = exp{—S) (localized states) and the solid line is the semi-Poisson distribution P{S) = ASexp{—2S) 
(critical states). The asymptotic distributions for extended and localized states are less relevant for our system since the data seem to 
approach the semi-Poisson law for the sizes used. 



The obtained spacing distribution shown in fig. 3 is neither Wigner (extended) nor Poisson (locaHzed). It is better 
described by a semi-Poisson intermediate level-statistics known to characterize pseudointegrable billiars Q , disordered 
systems at the critical point after performing certain average over boundary conditions [Q, quasiperiodic critical 
systems |8[^, interacting particles etc. Although the agreement with the semi-Poisson curve P{S) — 4:Sexp{—2S) 
is not perfect the obtained curve shows some scale-invariant characteristics and we can attribute the small disagreement 
with the semi-Poisson curve as due to the rather small number of levels considered. For example, the maximum size 
of 10946 considered in our computations, which gives the same number of levels, is clearly insufficient to answer 
definitely this issue, since over 100000 levels are usually required for this purpose [0,D, clearly impossible at present. 
However, we expect the agreement to improve as N increases. The discussion up to now also answers the second 
question. 
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FIG. 4. (a) The computed typical eigenvector amplitude distributions for the Fibonacci matrix of size = A^ig = 10946. (b) The 
corresponding scaling of the inverse participation ratio vs. size N for all eigenstate amplitudes, which defines their D2 fractal exponents 
to lie between 1/3 and 1. 



Finally, in order to answer the third question we have considered the eigenvector behavior of the Fibonacci matrix. 

Results for some typical eigenvector amplitude distributions together with the scaling of the inverse participation 
ratio are shown in fig. 4. We find that the vast majority of the eigenvectors are 'critical' sharing delocalized fractal 
characteristics. The scaling of the inverse participation ratio allows to compute D2 exponent for the eigenvector 
amplitudes. Wc obtain 1/3 < D2 < 1 while only very few eigenstates of this matrix are truly extended with D2 = 1 
usually lying at the spectral edges while states close to the band center are more 'localized' {D2 « 1/3). 

In summary, we have obtained numerically the full eigensolutions of a fractal Fibonacci matrix. The results show 
critical level-statistics with almost semi-Poisson nearest-level distribution and multifractal eigenvectors. The Fibonacci 
matrix studied shares similar characteristics to one-electron disordered systems at the mobility edge, quasicrystaline 
and also interacting electron systems. Moreover, it can be interesting for quantum systems which can show intermedi- 
ate level-statistics of what is believed to be 'critical chaos', in the border between quantum chaos and integrability. In 
might also offer some hints for the expected behavior from the diagonalization of matrices in complicated many-body 
Hamiltonians. 
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